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Abstract. We present two-dimensional simulations of wave propagation in a realistic, non- 
stationary model of the solar atmosphere. This model shows a granular velocity field and 
magnetic flux concentrations in the intergranular lanes similar to observed velocity and 
magnetic structures on the Sun and takes radiative transfer into account. 
We present three cases of magneto-acoustic wave propagation through the model atmo- 
sphere, where we focus on the interaction of different magneto-acoustic wave modes at the 
layer of similar sound and Alfven speeds, which we call the equipartition layer. At this layer 
acoustic and magnetic mode can exchange energy depending on the angle between the wave 
vector and the magnetic field vector. 

Our results show that above the equipartition layer and in all three cases the fast magnetic 
mode is refracted back into the solar atmosphere. Thus, the magnetic wave shows an evanes- 
cent behavior in the chromosphere. The acoustic mode, which travels along the magnetic 
field in the low plasma-/? regime, can be a direct consequence of an acoustic source within 
or outside the low-/3 regime, or it can result from conversion of the magnetic mode, possibly 
from several such conversions when the wave travels across a series of equipartition layers. 
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1. Introduction 

Classical helioseismology relies on waves with 
frequencies below the acoustic cut-off fre- 
quency. These type of waves are trapped inside 
the acoustic cavity of the Sun. High frequency 
waves above the acoustic cut-off frequency, 
however, are able to travel freely into the so- 
lar atmosphere. Along their path through the 
atmosphere these waves interact with the com- 
plex magnetic field that is present in the pho- 
tosphere and the chromosphere. Being able to 
observe these running waves in the atmosphere 
with high cadence instruments like MOTH 



dFinsterle et al.l l2004al) a new field in helio- 
seismology is emerging. Howev er, the analysis 
tools a t hand are still very basic dFinsterle et all 
l2004bt lHaberreiter etalj|2007l) and the char- 
acterization of the interaction of the magneto- 
acoustic waves with the magnetic field is still 
difficult. For the interpretation of the observa- 
tions, it is helpful to look at numerical sim- 
ulations of wave propagation through a mag- 
netically structured solar model atmosphere. 
In the past such simulations have always been 
carried out for idealized atm ospheric models 
and magnetic configurati ons dRosenthal et alj 
120021: iBogdan et alj|2003l) . However, in reality 
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the magnetic field is expected to have a com- 
plex structure. 

Here we report on two-dimensional simu- 
lations of magneto-acoustic wave propagation 
in a model of the solar atmosphere including 
realistic granular velocity fields and complex 
magnetic fields. In Sect. [2] we will give a short 
description of the setup of the atmospheric 
model and of the excitation of the waves. In 
Sect. [3] we will present and discuss snapshots 
of our wave propagation simulations. Results 
are summarized in the last section. 

2. Method 

For the simulations of magneto-acoustic wave 
propagation, we carry out two-dimensional 
numerical experiments with the CO s BOLD 
codefl The code solves the magnetohydrody- 
namic equations for a fully compressible gas 
including radiative transfer. Thus, the model 
atmosphere shows a realistic granular veloc- 
ity field and magnetic flux concentrations in 
the intergranular lanes caused by the advec- 
tive motion of the granules. More details on 
the application of the the CO s BOLD code for 
wave propag ation experiments can be found in 
ISteiner et all .20071) . 

2.1. Two-dimensional model of the solar 
atmosphere 

The model atmosphere used for the two- 
dimensional numerical wave propagation ex- 
pe riments is basically the same as was used 
bv lSteiner et al.l d2007l) . The computational do- 
main of the box covers the upper layer of the 
convection zone from about - 1 300 km below 
mean optical depth unity all the way to the 
middle layers of the chromosphere at about 
1 600 km above mean optical depth unity. 
Along the vertical direction an adaptive grid 
with 188 grid points is used corresponding to 
cell sizes of 46 km for the biggest cells in the 
convection zone and 7 km for the smallest cells 
in the chromosphere. Transmitting boundary 
conditions are applied for the lower and upper 

1 seewww.astro.uu.se/~bf/co5bold_main.html 
for the man pages of C05BOLD 
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Fig. 1. Magnetic setup of the model atmosphere 
represented in gray scales of log \B\. The white curve 
corresponds to the equipartition layer. On the right 
hand side a flux sheet with strengths up to 2 300 G 
has evolved. On the left hand side in the photosphere 
are weak horizontal magnetic fields. The crosses in- 
dicate the location of the wave excitation, where the 
numbers correspond to the three cases mentioned in 
the text. 



boundary. The lateral dimension of the box is 
5 000 km, with a grid resolution of 123 grid 
points corresponding to a horizontal spatial 
resolution of 40 km. Periodic boundary condi- 
tions are applied in the lateral directions. 

Figure Q] gives an impression of the mag- 
netic setup. The plot shows the absolute mag- 
netic field strength log |j5|. Since we are espe- 
cially interested in the interaction of the waves 
with the magnetic field, we use a model where 
a strong flux sheet has evolved during the sim- 
ulation time. In the core of the flux sheet, 
which is present on the right hand side of the 
box, magnetic field strengths up to 2 300 G 
are reached. On the left-hand side weaker hor- 
izontal fields have evolved. The white contour 
shows the equipartition layer where the Alfven 
speed equals the sound speed, and approxi- 
mately coincides with the /? = 1 layer. At this 
layer the acoustic and magnetic mode of the 
magneto-aco ustic waves can exchange energy 
(ICallvll2007l) . This layer is known as the mode 
conversion zone, where an acoustic wave can 
change to a magnetic mode and vice versa. 
The strength of this interaction depends on the 
angle between the wave vector and the mag- 
netic field (attack angle), the wave num ber k, 
and the width of the conversion zone (ICallvl 
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12007b . To investigate this interaction we excite 
waves at three different positions in the com- 
putational domain. The excitation positions are 
indicated by crosses in Fig. Q] The form of 
the excitation will be discussed in Sect. 12.21 
The location of the excitation is at a depth 
for which it was shown that the excitation of 
waves in the Sun is most likely to take place 
namel y between the surface a nd 500 km in 
depth dStein & Nordlundll200ll) . The three lat- 
eral positions are chosen at locations where in- 
teresting cases of wave and magnetic field in- 
teractions occur: 

i) where the magneto-acoustic waves interact 
with the magnetic canopy of the flux sheet, 

ii) where the magneto-acoustic waves are ex- 
cited inside the flux sheet, and 

iii) where the magneto-acoustic waves interact 
with weak horizontal magnetic fields. 




20 30 40 
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2.2. Wave excitation 

In order to excite waves in our simulation 
domain we follow the descripti on given by 
iParchevskv & Kosovichevl d2009l) . where the 
spatial and temporal behavior of the wave 
source is modeled by the function: 



f(r,t) = < 



where 



(1 -2T 2 )e~ T " for r<R src , 

for r > R src , 
(1) 



wo(f -to) • , 47T 
t = n with fo < t < to -\ , 

2 too 

R SIC is the radius of the source, A is the ampli- 
tude of the disturbance, and cog = 2nfo deter- 
mines the central frequency of the wave excita- 
tion. This function can now be used as a source 
function in the momentum equations by taking 
the gradient of Eq. Q] S = V/. Then, the com- 
ponents of the source function S can be com- 
bined with the components of the pressure gra- 
dient. Thus, the source function can be seen as 
a perturbation to the pressure gradient, launch- 
ing an acoustic wave with a central frequency 

Of (jjQ. 

For our simulations we use a central fre- 
quency of /o = 20 mHz and a radius of the 



Fig. 2. a) Temporal behavior of the wave source, b) 
Power of the wave source. The pressure disturbance 
excites a non-monochromatic wave with a central 
frequency of fo = 20 mHz. 



source of R src = 50 km. The temporal behav- 
ior of the wave source is plotted in panel a) of 
Fig. |2 The panel b) shows the power spectrum 
of the excitation source. A non-monochromatic 
wave is excited with most of the power around 
fo = 20 mHz. 

3. Results 

In this section we present snapshots taken from 
our simulations of wave propagation in a com- 
plex magnetically structured atmosphere. We 
will focus on the interaction between the wave 
and the magnetic field that is taking place in 
the mode conversion zone. 

Panel a) of Fig. [3] shows snapshots of the 
propagation of a wave that is excited just out- 
side of a flux sheet underneath the magnetic 
canopy. The propagation of both wave modes, 
the predominantly acoustic and predominantly 
magnetic mode, can be tracked in the pertur- 
bation of the square root of the energy den- 
sity 6e 1 ^ 2 = p l l 2 5v. Because of the position 
of the wave source beneath the (3 — 1 layer, 
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Fig. 3. a) Perturbation of the square root of the energy density, p i/2 6v, caused by a wave that is excited 
outside in the vicinity of the flux sheet. Shown are three snapshots taken at 80 s, 126 s, and 137 s into the 
simulation run. b) Perturbation SB of the magnetic field. The snapshots are taken at the same time as in 
panel a). In all plots, the solid contour corresponds to the equipartition layer indicating the conversion zone. 



in a region of high plasma B, the applied wave 
source excites a wave that is predominantly an 
acoustic mode. Upon reaching the conversion 
zone, indicated by the solid contour in Fig. [3j 
energy is transfe rred from th e acoustic to the 
magnetic mode (iCallvl 120071) . Because of the 
large attack angle between the local wave vec- 
tor and the magnetic field vector, most of the 
energy is converted to the fast magnetic mode. 
Transmission of the acoustic mode is not ap- 
parent. The magnetic mode is best seen in the 
perturbation 8B of the magnetic field, plotted 
in panel b) of Fig. [3] The second plot in panel 
b), at 126 s into the simulation, shows the fast 
magnetic mode emerging from the conversion 
zone. Because of the gradient of the Alfven 
speed, the wave is refracted back down into the 
photosphere. 

Next, we will focus on the wave propaga- 
tion when the wave source is located inside 



the flux sheet close to the equipartition layer. 
Figure |4] shows again snapshots of the pertur- 
bation of the square root of the energy den- 
sity 8e l/2 , panel a), and the perturbation 5B of 
the magnetic field, panel b). Since the wave 
source is located above the equipartition layer 
in the low-/? regime, the fast magnetic wave 
mode is excited directly and together with the 
slow acoustic mode, different from the previ- 
ous case. Because of the higher Alfven speed 
above the equipartition layer, the fast magnetic 
wave runs ahead of the slow acoustic mode. 
This fast magnetic mode can again be easily 
identified in the perturbation 6B of the mag- 
netic field, plotted in panel b) of Fig. [4] As in 
the previous case, the gradient of the Alfven 
speed causes the fast magnetic wave to re- 
fract to a degree that it travels back towards 
the lower atmosphere. The slow acoustic wave 
emerges from the location of the wave exci- 



4 



Nutto et al,: Numerical wave propagation in the solar atmosphere. 




IVrtiirfanrini A)B\ 



$d I 



sm iot» isoo ptnc e&pj xmo 3S«1 SOOS ■15O0 

[rrtO ■ BO S " |hn| 

tVrtitrbatlmi p'^Sv 



KM 1KG 1 500 HOC 2SOJ 3000 3500 4000 4500 
nm»» 105* *|ton] 

I Vi i ml i, in. in J iir 










f 



ism 20CQ 2500 0KO 4000 1500 

IViturlntiai in 



1 



1500 SOW 2900 3000 3500 4000 4900 



as 7T- 



DQQ 






&Q0 






b 






HQ 


r 




LQOC 




1 


500 1DO0 1500 2000 2SD0 30OTJ 3500 J 000 4500 


KW 


L ( t-ni|[l>arsi.ji f i ./i 





Fig. 4. a) Perturbation 6e [/1 caused by a wave that is excited inside of the flux sheet close to the equipartition 
layer. Shown are three snapshots taken at 80 s, 105 s, and 180 s into the simulation run. b) Perturbation SB 
of the magnetic field. The snapshots are taken at the same time as in panel a). In all plots, the solid contour 
corresponds to the equipartition layer indicating the wave-conversion zone. 



tation at a later time (last plot in panel a) of 
Fig.[4]i, and is then guided along the magnetic 
field lines into the higher layer of the atmo- 
sphere. As the acoustic wave travels up the 
atmosphere the wave front is steepening and 
starts to shock. 

Finally, we excite a wave just below weak 
horizontal magnetic fields. Snapshots of the 
wave propagation are plotted in Fig. [5] where 
panel a) shows again the perturbation of the 
square root of the energy density 6e 1 ^ 2 . Panel 
b) shows the perturbation of the pressure 
6plpo, with po being the local unperturbed 
pressure. As the excitation location is below 
the equipartition layer in the high-/? regime, 
primarily the fast acoustic wave mode is ex- 
cited. At the first equipartition layer, the lo- 
cal wave and magnetic vector are almost per- 
pendicular to each other. Thus, most of the 
acoustic wave energy is converted into the fast 



magnetic mode. A transmitted acoustic wave 
is not visible at first. However, the existence 
of multiple equipartition layers in the path of 
the magneto-acoustic waves, adds to the com- 
plexity of the interpretation. Multiple equipar- 
tition layers have never been addressed in for- 
mer investigations. Due to the large attack an- 
gle, one expects that no acoustic wave should 
be transmitted through the mode conversion 
zone. However, by looking at the pressure dis- 
turbance Sp/po in Fig. [5] it can be seen that 
an acoustic wave can be transmitted through 
the multiple equipartition layers by a series of 
mode conversion. In the middle plot of panel b) 
a strong pressure disturbance starts to reemerge 
above the multiple conversion zones. This can 
be seen as an indication of a transmitted acous- 
tic wave in spite of the horizontal magnetic 
fields. 
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Fig. 5. a) Perturbation 6e [/2 caused by a wave that is excited beneath weak horizontal magnetic fields. 
Shown are three snapshots taken at 80 s, 137 s, and 145 s into the simulation run. b) Perturbation of the 
pressure, Sp/po- The snapshots are taken at the same time as in panel a). In all plots, the solid contour 
corresponds to the equipartition layer indicating the wave-conversion zone. 



4. Summary 

We have shown numerical experiments 
of wave propagation in a realistic two- 
dimensional model atmosphere. We focused 
on the interaction between acoustic and mag- 
netic wave modes at the equipartition layer 
where the Alfven speed is similar the sound 
speed. 

When the location of the wave source is 
below the equipartition layer in the high-/? 
regime, the launched wave is predominantly 
acoustic in nature. In the first and third of the 
presented cases, the large attack angle of the 
local wave vector with the magnetic field vec- 
tor results in the conversion of the fast acoustic 
mode into the fast magnetic mode. When there 
are multiple equipartition layers, an acoustic 
mode can be transmitted into the higher layers 
of the atmosphere by a series of mode conver- 



sions. In all three cases the fast magnetic mode 
is usually refracted back towards the lower at- 
mosphere because of the gradient in Alfven 
speed. 

If the location of the wave source is close 
to the equipartition layer, both, the acoustic 
and the magnetic mode, are excited. The slow 
acoustic wave trails the fast magnetic wave. 
While the latter shows an upper turning point 
again, the slow acoustic wave is guided along 
the magnetic field into the higher layers of the 
atmosphere where the wave starts to shock. 

The simulations confirm that the propa- 
gation of magneto-acoustic waves in the so- 
lar atmosphere is strongly influenced by the 
equipartition layer. In comparison to investiga- 
tions where static background models are used, 
our simulations demonstrate that the equiparti- 
tion layer can show a complex topography and 
it is highly dynamic itself. 
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